A cryogenic surface-electrode elliptical ion trap for quantum simulation 
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Two-dimensional crystals of trapped ions are a promising system with which to implement quan- 
tum simulations of challenging problems such as spin frustration. Here, we present a design for a 
surface-electrode elliptical ion trap which produces a 2-D ion crystal and is amenable to microfabri- 
cation, which would enable higher simulated coupling rates, as well as interactions based on magnetic 
forces generated by on-chip currents. Working in an 11 K cryogenic environment, we experimentally 
verify to within 5 % a numerical model of the structure of ion crystals in the trap. We also explore 
the possibility of implementing quantum simulation using magnetic forces, and calculate J-coupling 
rates on the order of 10 3 s _1 for an ion crystal height of 10 /mi, using a current of 1 A. 



Quantum simulation is a technique whereby one may 
calculate properties of one quantum-mechanical system 
using a different, more controllable quantum system [1]. 
It has been shown, in principle, to hold great potential 
for solving problems in many-body quantum physics that 
are believed to be intractable on classical computers [2- 
5] . Recently, a particular interest has arisen in the simu- 
lation of quantum spin models using trapped ions [6] , and 
experiments on two and three ions confined in linear traps 
have been demonstrated [7, 8]. Direct analog quantum 
simulation of spin frustration could be performed using 
a two-dimensional (2-D) array of trapped ions. In ad- 
dition, such arrays have been proposed for more general 
quantum computations [9, 10]. 

Prior work toward developing 2-D arrays of ions has 
fallen within two paradigms: arrays of individual ion 
traps, and Coulomb crystals within a single trap region. 
Arrays of individual traps have been suggested as a way 
to implement simulations based on either optical forces 
[6] or magnetic gradient forces [11]. However, for re- 
alistic experimental parameters, both the motional and 
simulated interaction rates appear to be well below the 
expected decoherence rates in such ion traps [12], even 
when scaling to microscopic ion-ion distances. Despite re- 
cent advances in the design of arrays of surface-electrode 
ion traps [13], this scaling is likely to be problematic for 
some simulation protocols. Arrays of Coulomb crystals 
within a single trap, as realized in Refs. [14] and [15], 
bring other challenges, among them the E x B rotation 
of the ion crystal in Penning traps and rf heating due to 
micromotion in Paul traps. Recent theoretical work sug- 
gests that quantum simulation may be performed using 
2-D crystals of (9(100) ions in linear ion traps [16], even 
with the effects of rf and background gas heating. Al- 
though such traps may suffice for quantum simulation 
using optical forces, it is easier to achieve high mag- 
netic field gradients if the trap scale can be decreased. 
Magnetic interactions provide the additional advantages 
that errors due to spontaneous laser scatter are removed 
[11, 17] and that magnetic field gradients may be made 
uniform over the entire trap volume, enabling convenient 
global operations and mitigating errors due to micromo- 
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FIG. 1: Photograph of the elliptical ion trap. The rf elec- 
trode is labeled as Kf, while the others are dc electrodes. 
The rf electrode is approximately elliptical; however, in order 
to tilt the trap principal axes and improve Doppler cooling, 
the semiminor axis of the electrode is stretched along one di- 
rection. The semimajor axis is labeled B in the figure, while 
the semiminor axes are labeled A and A 7 . The semimajor 
and semiminor axes of the center electrode are 0.94 mm and 
0.71 mm, respectively. The coordinate system is defined in 
the figure; the z axis points out of the page. 

In this Letter, we present the design and proof-of- 
principle demonstration of a Paul trap that generates a 
2-D crystal of trapped ions and that is amenable to mi- 
crofabrication. Our trap is a surface-electrode ring trap, 
similar to the trap presented in Ref. [18], but with the 
usual cylindrical symmetry broken in two ways. The first 
is that the center ground and rf ring electrodes are based 
on ellipses rather than on circles, similar to the ellipti- 
cal trap invented by DeVoe [19], in order to break the 
degeneracy of the in-plane motional modes and give the 
ion crystal a unique orientation in space. The second 
is that the rf ring electrode is made wider on one side 
than the other, to tilt the ion crystal with respect to the 
plane of the ion trap and improve Doppler cooling. We 
first present a numerical model of the trap, along with 
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experimental measurements which verify the model. We 
then turn to the question of what conditions must be met 
for a quantum simulation to be performed within such a 
trap. To that end, we explore theoretically how quantum 
simulations may be done using magnetic field gradients 
in a scaled-down version of the elliptical trap. 

The trap is a printed circuit board consisting of cop- 
per electrodes on a low-loss Rogers 4350B substrate, the 
same manufacturing process that was used in Refs. [20] 
and [21]. A photograph of the trap is presented in Fig. 1. 
The trap was numerically modeled using the Charged 
Particle Optics (CPO) boundary-element software pack- 
age. Typical secular frequencies for an rf drive frequency 
of ft/(27r) = 3.5 MHz and an rf amplitude of V rf = 150 V 
are v x = 124 kHz, v y = 110 kHz, and v z = 234 kHz. 
Confinement may achieved by applying an rf voltage and 
grounding all dc electrodes. 

A cryogenic environment was used for our experiments, 
which offers the advantages of suppressed electric field 
noise [22] and a rapid pumpdown to ultra-high vacuum 
pressures. The trap was mounted in a closed-cycle cryo- 
stat (ARS model DE-210 with DMX-20 anti-vibration 
interface), similar to one described in Ref. [23]. The elec- 
trode labeled V3 was thermally anchored to the base- 
plate, which typically reached a temperature of 11 K. 
The vacuum pressure was upper-bounded by reading the 
current on a 20 1/s ion pump that was connected to the 
room-temperature vacuum housing. The vacuum pres- 
sure at room temperature was at most 5 x 10 -10 torr, 
and within the cryogenic environment was expected to 
be much lower. 88 Sr + ions were loaded by photoioniza- 
tion of a neutral strontium beam and Doppler cooled on 
the 422 nm 5Si/2 —> 5Pi/2 transition with a 1092 nm 
repumper. Typical powers were 50 /iW of 422 nm and 
100 /iW of 1092 nm with beam waists of « 50 /im. 

For the measurements presented here, an rf amplitude 
of Vrf = 150 V at Q/(2tt) = 3.5 MHz was used. Stray 
fields were compensated by adjusting the dc voltages such 
that a cloud of ions remains stationary when V r / is varied 
[24]. The final set of voltages used was V\ = —3.90 V, 
V 2 = 1.56 V, V 3 = V A = V, and V c = -2.62 V. Sec- 
ular frequencies were then measured by applying a low- 
amplitude resonant signal to one of the dc electrodes, 
and recording the dip in the ion cloud's fluorescence as 
the ions were resonantly heated. For the above condi- 
tions, we measure v x = 177 kHz, v y = 141 kHz, and 
v z = 414 kHz, with an error « 2 kHz on each. These 
frequencies differ from those listed above due to the fact 
that the dc voltages used for compensation also alter the 
curvature of the trapping potential. 

The structure of crystals in the trap is calculated by 
numerically minimizing the potential energy of a set of TV 
ions in a harmonic potential characterized by the three 
secular frequencies. For sufficiently small iV, the crystal 
extends only along x and y. For v x « v y , the condition 
for a planar crystal may be written as v z > y/70Ni^ /ir 3 , 
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FIG. 2: Images of crystals of two and four ions in the elliptical 
trap. The ion-ion spacings d y and d x are along y and x, re- 
spectively. Since the separation of two ions may be calculated 
analytically, the two-ion crystal is used to calibrate the mag- 
nification of the optics system. For two ions, d y = 16.5 /mi. 
For four ions, d y = 28 =b 3 /im and d x = 17 =b 3 /im, compared 
to the calculated values of d y = 29.6 /im and d x = 17.1 /im. 
Note that d x 7^ d y due to the eccentricity of the trap. 



as noted in Ref. [25]. For v x ^ v y , such a condition 
may be computed numerically. Fig. 2 contains images 
of crystals of two and four ions, along with their mea- 
sured separations. We find that for the measured secu- 
lar frequencies, the agreement is quite good. Thus, our 
model of the elliptical trap potential and crystal structure 
therein is validated, at least for small numbers of ions. A 
limited 422 nm laser stability, along with the possibility 
of rf heating, are possible technical explanations for why 
larger ion crystals were not obtained. Larger crystals 
were later observed in a similar ion trap (with circular 
geometry) in the same cryogenic apparatus [18]. 

This trap design provides a unique platform for quan- 
tum simulation, in that wires that produce magnetic field 
gradients may be integrated into or beneath the ground 
electrode (V c ). Although quantum operations based on 
oscillating magnetic fields have been proposed [26], we 
concentrate here on static field gradients. For a given 
configuration of current- carrying wires, the magnetic 
field B above the trap is calculated by direct numerical 
integration of the Biot-Savart law. The state-dependent 

force F on a trapped ion is given by F = — V ^/2 • B^j , 

where jl is the magnetic moment of the ion. A qubit may 
be formed from either a Zeeman-split electronic ground 
state or a pair of hyperfine levels. 

As an example of how a quantum simulation may be 
implemented using this system, we calculate the mag- 
netic fields for three concentric square wires embedded 
within the center ellipse. Each is assumed to carry a 
current of 1 A, on par with the currents used in neutral 
atom chip traps [27]. The wire configuration is depicted 
in Fig. 3. Using these calculated force values, and as- 
suming the orientation of the atomic magnetic moment 
is along y (as set by an external bias field), one may ob- 
tain an order of magnitude estimate for the simulated 
interaction rate of a quantum simulation. We choose to 
calculate the effective J-coupling rate for a simulation 
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FIG. 3: Left: Example wire configuration in the plane of an 
elliptical ion trap. Right: A table of some values of the simu- 
lated coupling rate J as a function of the ion height H above 
the trap surface, for 88 Sr + . The ion height is varied by scaling 
down the entire trap structure, including the current-carrying 
wires, and the corresponding approximate ion-ion spacing d 
is given for reference. We assume a current of 1 A through 
each wire. The secular frequencies are scaled up in inverse 
proportion to the trap scale. The quadratic dependence of J 
on F indicates that these values of J could be improved by a 
modest increase in the current or in the number of concentric 
wire loops. 



of quantum spin models; according to Ref. [6], and for 
two ions undergoing an antiferromagnetic simulated cou- 
pling, this rate is given by hJ = Ke 2 c F 2 / (647r 5 eora 2 ^ 4 <i 3 ), 
where k is a constant of order unity, F is the magnitude 
of the state-dependent force, m is the ion mass, v is the 
relevant secular frequency, eo is the permittivity of free 
space, e c is the ion charge, and d is the mean distance 
between neighboring ions. Some calculated values of F 
and J as a function of the overall trap scale are given in 
Fig. 3, showing that a trap scale on the order of tens of 
microns can produce (9(10 3 s _1 ) coupling rates. 

Although sufficient interaction rates for the simulation 
of quantum spin models are theoretically possible in our 
trap, there are some other effects that may reduce the 
fidelity of a quantum simulation that must be addressed. 
These include heating of the ion crystal, micromotion, 
and variations in the distance between ions. The heating 
rates of planar ion crystals in linear ion traps have been 
calculated in Ref. [25] ; the authors predict that for up to 
(9(100) ions, high-fidelity quantum simulations should be 
possible despite modest heating. In a 2-D crystal, some 
micromotion will be present, but micromotion represents 
a breathing-mode oscillation, in which the relative micro- 
motion amplitude between nearest-neighbors is much less 
than the total amplitude of micromotion. Furthermore, 
magnetic gradient-based forces vary much less rapidly in 
space than optical forces, which mitigates the effects of 
micromotion. Finally, we note that the ion-ion spacing is 
quite uniform near the center of a many-ion crystal, pro- 
vided v x « v y . Further calculations indicate that, with 
appropriate dc and rf voltages, a planar crystal of (9(100) 
ions could be produced in our trap [28] . Assuming perfect 
quantum operations, this number would in principle be 
enough to outperform exact classical simulations. While 



the possible number of ions in a planar crystal is com- 
parable to the state of the art before this work, our trap 
design offers the possibility of harnessing the benefits of 
microfabricated ion trap technology, including magnetic 
forces, for quantum simulation. 
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